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ABSTRACT 



FUSE, HST and SDSS spectra of the cataclysmic variable SDSSJ080908.39 
+381406.2 provide a spectral flux distribution from 900-9200A. This data set 
is used to illustrate procedures for calculating and testing system models. The 
spectra are not contemporaneous. The illustrations are based on a system with 
a I.OMq white dwarf, a O.SOM©, 3500K, Roche lobe-filling secondary star, and 



^ Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc. under NASA contract NAS5-26555, the NASA-CNES-CSA Far Ultraviolet Explorer, which is operated 
for NASA by the Johns Hopkins University under NASA contract NAS5-32985, and with the Apache Point 
Observatory 3.5 m telescope, which is operated by the Astrophysical Research Corporation. 
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an accretion disk extending to the tidal cutoff radius. Assuming a similar accre- 
tion state for the non-simultaneous spectra, the best standard model fit is with a 
mass transfer rate of 3.OxlO~^M0yr~^. Extensive simulations demonstrate that 
the accretion disk must be truncated at its inner edge if the temperature profile 
follows the standard model, but truncated models face severe objections, which 
we address. Following additional simulation tests, we obtain a model accretion 
disk with a temperature profile comparable to the profile for SW Sex as deter- 
mined from tomographic image reconstruction. This model fits the discovery 
SDSS spectrum well but has a flux deficit in the UV and FUV. Emission from a 
white dwarf is a plausable source of additional flux. Adding this source to the 
disk synthetic spectrum produces FUV flux that can explain the observed flux. 

An additional (archival) SDSS spectrum is fainter by about 0.3 magnitude in 
the optical. Additional analysis showed that UV residuals from a model fitting 
the archival optical wavelength spectrum are unacceptably large. Contemporane- 
ous spectra from all wavelength regions would be necessary for a reliable system 
model. Our discussion illustrates how this conclusion follows from the system 
models. 

Subject headings: accretion, accretion disks — novae, cataclysmic variables — 
stars: individual(SDSSJ0809) — ultraviolet: stars — white dwarfs 



Introduction 



Cataclysmic variables (CVs) are semi-detached binary stars in which a late main se- 
quence star loses mass onto a white dwarf (WD) via Roche lobe overflow. In systems con- 
taining a non-magnetic WD, accretion proceeds through a viscous disk. However, in systems 
containing a magnetic WD, accretion may proceed directly from the inner Lagrangian point 
onto field lines in the strong field case, or through a partial disk in the case of interme- 
diate field strength. The group of CVs that have typical orbital periods between 3-4 hrs, 
strong single-peaked Balmer emission, strong Hell emission and deep central absorption 
in the Balmer and Hel lines (usually n e ar orb i tal phase 0.5) ar e calle d the SW Sex stars 



(e.g.. iHoneycutt. Schlegel. &: Kaitchuckl (119861 ): ISzkody &: Pichg (Il990l ): iThorstensen et al. 



(Il99ll )). As a group, their orbital periods place the SW Sex stars at the top of th e period 
gap, which likely is an aspect of the normal secular evolution of the CV population (jWarner 
19951 ). The SW Sex stars typically have high mass accretion rates and a f ew of them may have 



interm edia te streng t h ma gnetic fields as evinced from weak polarization (IRodriguez-Gil et al. 



200ll ). See IWarnerl (119951 ) for a thorough review of CV types and behaviors. 
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The observational data that are commonly available to study CVs include spectra and 
photometry. Deduction of the physical properties of these binary stars depends on simulation 
of the observations, based on a physical model believed to represent reality with adequate 
accuracy. The accretion disk is typically the most important radiating component in a CV 
over a wide range of wavelengths. A comrn o n practice h as been to use synth etic stellar spectra 
to sim ulate the accretion disk (lla Dous I (Il987l . Il994l )). More recently, IWade fc Hubeny 
( 1l998l ) have calculated a set of models based on the TLUSTY program. It is desirable 
to have a simulation program that calculates an accretion disk ab initio, given an orbital 
i nclination, W D mass, and m ass transfer rate. A com bination of the programs TLUSTY 
(lHubenylll988l ) and BINSYN (ILinnell fc Hubenylll996l ) provides this feature. The details of 
how the calculation takes place, in the context of a specific system, have not been described 
in the literature. The availability of both HST and FU SE data combined with a good SDSS 
spectrum of a bright system (few systems provide this combination) provides the impetus 
for modeling this system. Since the data were not contemporaneous, a definitive solution is 
not available. However, as we show later, we are able to find a solution which satisfactorily 
represents the data set. At the same time, the data set provides a very useful basis to 
illustrate the details of building and testing system models, and the results turn out to be 
of considerable intrinsic interest. 



2. The SDSSJ0809 system 

The Sloan Digital Sky Survey (York et al. 2000) provides 5 color photometry (ugriz) and 
spectra (3800-9200A at a resolution of 3A) of the survey objects. The survey has identified 
a large number of CV systems (Szkody et al. 2002, 2003, 2004, 2005). 

During the survey, SDSS J080908.39+381406.2 (hereafter abbreviated as SDSSJ0809) 
was discovered as a bright CV (g~15.6) (Szkody et al. 2003). The discovery paper provides 
the photometric colors and optical spectrum of SDSSJ0809 as well as the results of time- 
resolved followup spectra obtained over the course of 3 hrs. The SDSS spectrum shows a 
steep blue continuum with strong Balmer emission lines as well as prominent emission from 
He II A4686 and the CN blend at 4640A. The time-resolved spectra show deep and transient 
absorption in the Balmer and Hel lines. If the Ho; and H/3 radial velocities indicate orbital 
phases, then the deep absorptions occur near orbital phase 0.5 (Szkody et al. 2003), which 
is a hallmark of the SW Sex stars. However, we lack any confirming evidence (eclipses) and 
we are presently unable to make any definite statement about orbital phases. There was no 
noticeable eclipse during the 3 hrs of observations and the radial velocity curve constructed 
from the Ha and H/3 emission lines indicated an orbital period near 2.4 hrs. However, the 
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shortness of the data interval precluded the determination of an accurate period. Recent 
photometry over many nights indicates that the period is 192 min (more typical for SW Sex 
stars) and there may be a shallow or partial eclipse (Rodriguez- Gill et al., in preparation) 
with amplitude of less than about 0.1 magnitude. We will use the 192 min period throughout 
this paper. 

Due to the brightness of this system, ultraviolet measurements were possible with both 
the Far Ultraviolet Spectroscopic Explorer (FUSE) and the Hubble Space Telescope (HST) 
and were used to explore the accretion parameters and geometry of the high accretion regime. 
In the high inclination SW Sex star DW UMa, which shows deep eclipses, the accretion disk 
completely obscures the WD and inner disk during high states of accretion (Knigge et al. 
2004). When the accretion is reduced during l ow states of mass transfer, th e underlying WD 



is revealed and it is very hot (near 50,000K (lAraujo-Betancor et al.ll2003l )). Most SW Sex 



stars have high inchnations and deep eclipses such that, as in DW UMa, the WD is hidden 
from view by the accretion disk. We observed SDSSJ0809 in the hope that its apparent 
lower inclination (as evinced by the lack of deep eclipses) might permit a more direct view of 
the WD and inner disk, and thus allow us to explore the characteristics of these components 
of an SW Sex star. An objective of this investigation is to determine the extent to which 
the large SDSS spectral range, combined with spectra from HST and FUSE, together with 
related collateral information, permit development of a self-consistent model. 

As the following sections indicate, we have a single HST spectrum, and two sets of SDSS 
and FUSE observations. None of the observations are simultaneous. The photometric ob- 
servations Rodriguez- Gill et al. were obtained some 300 days after the HST exposure and 
continued over an interval of 26 days. The photometry shows only minor orbital variation, 
with an amplitude of 0.1 magnitude that may indicate a grazing eclipse, but with no sys- 
tematic trend. There is consistent overlap of the HST and one of the FUSE observations. 
Since the SDSS photometry, the first SDSS spectrum, the subsequent APO spectroscopy, the 
Rodriguez-Gil photometry and the optical measurement from the HST snapshot all agree 
(to ±0.1 mag), we first treated the system as if it is constant at all observed wavelengths 
longward of the FUSE observations. 



3. FUSE Observations and Data Processing 

We observed SDSSJ0809 with FUSE during two time intervals approximately 1 year 
apart; Observation 1 occurred on 28 Mar 2003 and Observation 2 occurred on 16-17 Mar 
2004. Because of the FUSE orbital constraints, each observation was divided into several 
exposures (4 exposures for Observation 1, 8 exposures for Observation 2) - see Tabled] for 
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the FUSE observing log. All data were obtained using the LWRS aperture an d TTAG accu- 



mulat ion mode (for FUSE spacecraft and instrument details see, for example, ISahnow et al. 
200dl^ I). We used the CalFuse v3.0.7 pipeline software to prepare intermediate data files from 
the raw data files obtained during each exposure. We then used the IDL routine cf_edit 

v2.sE, as well as custom-built IDL routines, to extract spectra from the various detector and 
mirror segments and combine them into a time-averaged spectrum for each observation. The 
time-averaged spectra have total equivalent exposure times of ~ 8.4 ksec for Observation 1 
and ~ 10 ksec for Observation 2 (total exposure times from the mirror /detector pairs differ 
by up to ~ 5% due to differences in rejected data; the combined spectrum accounts for this 
by weighting by exposure time when combining data from different mirror/detector pairs). 
The final combined spectra (see Figures [1] and [2]) were rebinned onto a uniform wavelength 
scale with dispersion 0.20 A pixel"^ by averaging flux points from the original dispersion 
(0.013 A pixel"^) into wavelength bins of width 0.20 A. 



4. Analysis of FUV Spectra 

Our initial plan was to combine the two FUSE data sets to improve the S/N of the final 
time-averaged spectrum. However, it was apparent that the spectra obtained from the two 
observations were somewhat different, so we treated them separately. Shortward of lOOOA 
the two spectra are approximately equivalent, although this apparent agreement is likely 
due mainly to the weak detection at these short wavelengths. None of the sul fur or silicon 



emission lines often seen i r i FUV spectra of CVs in this wavelength region (e.g., iHoard et al. 



2003; 


Fronine et al. 


2001; 


Long et al. 


1994) 



of Figures [T] and [2]) . The C III, N III emission complex is possibly present in both spectra; 
however, it coincides with the wavelengths of a series of airglow features, which makes a 
definitive detection difficult. 

At longer wavelengths, the Observation 2 spectrum is approximately twice as bright in 
the continuum as the Observation 1 spectrum; for example, the "mean" continuum level at 
1155A is ~ 0.75 X 10~^^ erg s~^ cm~^ A~^ in the Observation 1 spectrum and 1.5 x 10"^'* 
erg s~^ cm~^ A~^ in the Observation 2 spectrum. 

There are also differences in the ionic features present at A > lOOOA in each spectrum. 
In the Observation 1 spectrum, the S-I-Si line complexes at 1055-1085A and 1105-1130A are 



^Also see the FUSE Science Center web page at http : //fuse . pha ■ jhu . edu/[ 

^Written by Don Lindler and available at http://fuse.pha.jhu.edU/analysis/fuseJ.d.l_tools.html 
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present as broad, undifferentiated "humps" in the continuum. However, in the Observation 
2 spectrum, these same wavelength regions contain clearly differentiated emission and/or 
absorption features. The most striking difference is in the C+N+0 line complex at 1165- 
1185A. In the Observation 1 spectrum, this is present mainly as a fairly strong emission 
feature (strongest of the non-airglow lines) at the wavelengths spanned by the O III multi- 
plet. In the Observation 2 spectrum, on the other hand, this line complex displays a profile 
consisting of a broad hump from 1165-1185A, with a possible discrete O III emission compo- 
nent, and deep absorption (reaching down to at least the level of the underlying continuum) 
at the wavelengths spanned by the C III+N I multiplets. 

The O VI lines at AA1031.9, 1037.6A, which are often among the mo st prominent emis- 



sion features in the FUV spectra of CVs with ongoing accretion (e.g., iHoard et al.l 12005 



20081 . 12OO2I : lMauchel[l999l ). are only weakly present in both spectra of SDSSJ0809. In addi- 



tion, the line profiles of O VI appear to be in emission on the blue side and absorbed on the 
red side (similar to the shape of the C+N-l-0 line complex at 1165-1185A in Observation 2). 
Detection of the S III multiplet at 1012-102lA is even weaker (or non-existent) compared to 
O VI in both spectra. Neither spectrum displays detectable Si III features at 1144A. 

The total exposure time for each of the time-averaged FUV spectra of SDSSJ0809 is 
shorter than the orbital period of the CV (75-85% of Porb = 192 min). However, the 
total orbital phase coverage during a given observation differs from the other observation by 
only about 10% and spans most of an orbit. Thus, we beheve that the differences between 
the spectra are not likely to be due only to differences in lines-of-sight through the orbital 
geometry. In any case, the current estimate of the orbital period of SDSSJ0809 is not known 
precisely enough to preserve cycle count well enough from Observation 1 to Observation 2 
that we could compare relative phases of the two FUSE spectra. 

The Observation 1 spectrum, which we designate as FUSEl, accords with the HST 
spectrum in the wavelength overlap region better than the Observation 2 spectrum. The 
Observation 2 spectrum difference from the Observation 1 spectrum could arise from a 
change in the accretion disk (discussed subsequently) in the time interval between the two 
observations. 



5. HST Observations 



SDSSJ0809 was observed on 2004 Apr 29 as part of aii HST Spa.ce Telescope Imaging 
Spectrograph (STIS) snapshot program (see iGansicke et al.l (120031 ) for a description of this 
program). Using the STIS/CCD acquisition image, we determined that the magnitude in 
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the F28x50LP filter (~ 5500 - 11,000A) was about 15.5. The STIS spectrum is an 800- 
sec exposure obtained with the G140L grating, in accum mode, which provided wavelength 
coverage from 1100-1700A with a resolution of 1.2A. The data were reduced with CALSTIS 
(V2.13b) which corrected for the decaying sensitivity of the grating. 

The spectrum is shown in Figure 3. The continuum displays a downturn near Lya and 
weak emission lines of GUI, NV, CII, SilV and CIV can be identified. The lines are noticeably 
weaker than those in the high state STIS spectrum of DW UMa (Knigge et al. 2004), 
consistent with a lower inclination for SDSSJ0809. The absorption cores at the SilV doublet 
wavelengths, superposed on the overlapping doublet emission, and the doubled appearance 
of the CIV lines are similar to the STIS spectrum of DW UMa at phase 0.5 (Knigge et 
al. 2004), when the absorption events characteristic of SW Sex stars are most prominent. 
A number of additional line identifications in Figure 3 mark absorption cores associated 
with broad emission features; in most cases the absorption core goes below the continuum 
level, indicating a substantial column density along the line of sight. Unfortunately, we 
cannot determine the phase of the SDSSJ0809 STIS spectrum as there currently is no orbital 
ephemeris. 



6. SDSS and APO Observations 



The discovery spectrum of SDSS0809 from MJD52251 is described in ISzkody et al. 



(120031 ) and is reproduced here in Figure 4. The exposure time, 5803 sec, spans approximately 
half an orbital cycle. There is no indication of a Balmer discontinuity. The SDSS photometry 
for SDSSJ0809 gave g=15.61, and the 5500A flux in the SDSS spectrum gave an approximate 
V magnitude of ~ 15.4. 

The Apache Point Observatory (APO) follow- up spectra, obtained over 2.9 hours, pro- 
vided a radial velocity curve in both Ha and Rp JSzkodv et al.ll2003h . The time-resolved 
APO spectra showed the strong modulation and deep transient absorption in the He I and 
Balmer line cores that characterize SW Sex systems ( Szkody et al.lEo03l ). He II and CN 
emission are present at all phases of the APO spectra. 



7. CV accretion disks 



Frank. King. &: Raind (119921 ) discuss the theory of accretion disks. They show that the 
effective temperature of an accretion disk, in hydrostatic equilibrium and timewise invariant, 
at a distance R from the accretion disk center, is given by the equation 



- 8 - 



Tefr(i?) 



(3GMM 
I 87ri?V 



1 - — 



R, 
R 



1/2 



1/4 



where M is the mass of the WD, M is the mass transfer rate, and i?* is the radius of the 
WD. G is the gravitation constant and a is the Stefan-Boltzmann constant. This equation 
defines the standard model for the radial temperature profile of a CV accretion disk. 

The accretion disk is heated by viscous shear. It is now believed that the p hysical 
basis for viscosity is the magnetic instability mechanism (IBalbus fc Hawley 1119911 ). Lack 
of knowledge of a physical basis for viscosity in the original study of accr etion disks led 
to introduction of the dimensionless parameter a (jShakura fc Sunyaevlll973l ). An alternate 
viscosity presc ription is to use a Reynolds number, typically indicated by the symbol Re. 
Pringld (Il98ll ) discusses both prescriptions and shows that, for CV accretion disks that 
undergo outbursts to have temporal characteristics roughly equal to observed values, then 
the a prescription requires values in the range 0.1-1.0, and the Reynolds number prescription 
should use values in the approximate range 10^-10'^. Although important a dvances have been 



made i n specifying the vi s cosity profile within an a c cretion d i sk an nulus (IStone fc Balbus 



fll996f ): IStone et all (Il996f ): iBalbus fc Hawlev I (119971 ): iBalbus I (j2002|)), it still is not possible 
to calculate a value of a or Re from fundamental physical principles to use in modeling a 
particular annulus. An important point is that equation (1) has no explicit dependence on 
the viscosity parameter. We assume that a single value of the viscosity parameter applies to 
an entire accretion disk. 

We mod el a CV a ccreti on disk as a nested series of cylindrical annuli. The theoreti- 



cal model by iHubenyl (jl990l . he reafter H90) provides the model for the individual annuli 
with additional 



Kfi'z & Hubenv 


(11986). Computer 


( Hubenv 


1988; 


Hubenv & Lanz 


1995) 



TLUSTY models uses program SYNSPEC (v. 48) (IHubeny. Stefi fc Harmanedll985l : iHubeny 



19901 ). An annulus synthetic line spectrum includes physical line-broadening effects, such 
as natural, van der Waals, Stark, or from turbulence. We adopted solar composition for all 
synthetic spectra. A useful option in SYNSPEC is to calculate continuum spectra, but also 
including H and He lines, rather than line spectra. 

H90 develops the theory of the vertical structure of an annulus. The vertically-averaged 
viscosity and its specific value as a function of reference level within the annulus appear 
explicitly in this theory. The annulus Teg as a function of radial distance from the geometric 
center is given by equation (1). 



Calculation of a model for a single annulus proceeds in four steps: (1) calculation of 
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a gray model using mean opacities rather than frequency-dependent ones; (2) iteration, 
starting with the gray model but now using frequency- dependent opacities, while imposing 
local thermodynamic equilibrium (LTE) excitation relations, until changes in a state vector, 
at all levels, become less than a prescribed small value, producing a converged model in 
LTE; (3) relaxing the LTE requirement and iterating to produce a non-LTE model for the 
continuum; (4) using the model from (3) and including line opacities and again iterating to 
produce a final non-LTE model. 

Our experience is that models for inner annuli typically converge but convergence be- 
comes more and more fragile for larger radius annuli until only a gray model can be calcu- 
lated. Depending on the mass transfer rate, even the gray model calculations may fail for 
the largest radius annuli. Based on extensive tests, we find that plots of annulus synthetic 
spectra progress smoothly from smaller to larger annulus radii, even though the models 
transition from LTE to gray models. Consequently, to minimize the required amount of 
computer time, we have used gray models exclusively in this investigation. 



7.1. Calculation of synthetic spectra for a CV system 



The BINSYN software suite ( iLinnell fc Hubenyi |l996|) is used to calculate synthetic 



spectra for comparison with the observational data. For an accretion disk system, BINSYN 
produces separate output synthetic spectra for the total system, the mass gainer, the mass 
loser, the accretion disk face, and the accretion disk rim (including a hot spot, if specified). 
Line-broadening and line displacement due to the Doppler effect are calculated in BINSYN 
as the programs generate synthetic spectra for the total system and the individual system 
components. The calculated synthetic spectra require prior assembly of arrays of source 
synthetic spectra produced externally to BINSYN, one array for each object (mass gainer, 
accretion disk face, etc.). Thus, for a model corresponding to Table 2 (see below) and with 
a 3500K secondary star, BINSYN would require 21 synthetic spectra corresponding to the 
annulus radii in Table 2, a synthetic spectrum for the WD, a synthetic spectrum (in this 
case duplicating the outer annulus spectrum) for the accretion disk rim, and a synthetic 
spectrum for the secondary star. For the latter object, since its contribution in the present 
case is negligible, a single synthetic spectrum meets the formal requirement. Otherwise, 
because the secondary star fills its Roche lobe, a range of synthetic spectra would be required 
covering the photospheric variation in T^s and log g. The BINSYN model for the distorted 
secondary star includes effects of gravity darkening and irradiation by the primary and the 
accretion disk. Shadowing of the secondary, by the accretion disk, from irradiation by the 
primary is included. 
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The number of annuli specified, e.g., in Table 2, must be large enough that interpo- 
lation to a non-tabular T^s value can be accurate. The output spectrum is produced by 
interpolation within the array of synthetic spectra to produce, in effect, a local photospheric 
model atmosphere at each photospheric segment. Integration over the segments, with proper 
allowance for Doppler shifts and sources of line broadening, and with suppression of con- 
tributions from segments hidden from the observer, then produces the synthetic spectrum 
for the object (star, accretion disk, rim) at the particular orbital inclination and longitude 
under consideration. The system synthetic spectrum is the sum of contributions from the 
separate objects. Stellar objects are represented by the Roche model, including allowance 
for rotational distortion up to critical rotation. 

The accretion disk is represented within BINSYN by a specified numb er (typically 32) of 



concentric annuli with fixed radial width but of increasing thickness (see iLinnell &: Hubeny 



(Il996l )). up to a specified rim height. Each annulus is divided into azimuthal segments (typ- 
ically 90). Note that the number of annuli specified within BINSYN usually exceeds the 
number of calculated annulus models. This feature is necessary to provide adequate resolu- 
tion in calculating, e.g., eclipse effects. BINSYN calculates a Tes va lue for each internally 



speci fied annulus (the 32 mentioned above); these are standard model (IFrank. King, fc Raine 



19921 ) values by default, but an option permits assigning an individual T^s to each annulus. 
This feature allows evaluation of an arbitrary temperature profile for the accretion disk, 
including an isothermal model. The program optionally permits calculation of irradiative 
heating of the accretion disk by the WD, based on a bolometric albedo formalism. The 
output spectra from BINSYN are in the same mode as the annulus, etc., spectra-either 
line spectra or continuum spectra (i.e., if the synthetic spectra for the annuli are continuum 
spectra, that mode having been chosen for the SYNSPEC calculations, then the system, etc., 
synthetic spectra produced by BINSTN also will be continuum spectra). 

BINSYN requires input specification of the inner and outer radii of the accretion disk. 
This feature provides important flexibility to truncate the disk at an inner radius and to 
set the outer radius at the tidal cutoff radius as dictated by the mass ratio of the stellar 
components. The accretion disk temperature profile typically will place all of the accretion 
disk on the "ho t branch" re quired for dwarf novae outbursts by the Disk Instability Model, 



hereafter DIM flOsakilll996f ) 



8. WD and secondary star masses and T^s values 



The observed spectra of SDSSJ0809 give no immediate indications of an underlying 
WD spectrum, apparently preventing a direct determination of the WD Te^. The sub- 
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sequent discussion shows that SDSSJ0809 has an accretion disk, so it is clear that the 
secondary star m ust fill its Roche lobe. R e cent y ears have seen maj or progress in rnodels 



of low-mass stars (IBaraffe et al.l Il995l . 119971 . Il998l ) and recent papers (IKolb &: Baraffd 12000 



Andronov &: Pinsonneaultll2004l ) detail the applica t ion of these models to the secondary stars 
in CVs. From Figures 2 and 3 of iKolb fc Baraffd (120001 ) we determine a secondary compo- 
nent spectral type of about M3 and a secondary component mass M2/M0 = 0.30 ± 0.10. 
This determination has used the theoretical ZAMS, and there is a substantial difference 
between observationally determined spectral types and the theoretical ZAMS in the period 
range 2 — S'^. This difference leads to the large uncertainty in the secondary component mass. 
Subsequent sections consider the propagation of this uncertainty into the system geometry. 



Smith fc DhillonI (119981 ) specify an average WD mass of 0.69Mq for WDs below the 
period gap and O.SOMq for WDs above the period gap. Their tabulation includes only one 
object in the period interval 2^ to 4^, and it is a polar. The systems with periods between 
4^ and 5*^ have WD masses between O.GM© and I.26M0, with unreliable values in several 
cases. The 192 min period places the system above the period gap. For the purposes of an 
initial illustration we adopt a WD mass of I.OMq, and with an uncertainty of O.2M0. 

The lack of deep eclipses sets an upper limit on the orbit al inclination. The p resence of 
absorption cores associated with emission features (Figure 3).( Szkody et al.lbood ). with the 



absorption cores going below the continuum level, sets a rough lower limit on the inclination. 
We have determined the orbital inclination for which the secondary component barely fails 
to eclipse the accretion disk rim, both for the q = 0.30 case and for a g = 0.20 case. For 
the q = 0.30 case, the inclination is i = 65°, and for the q = 0.20 case the inclination is 
i = 67°. We take the observed low-amplitude light curve, together with the absorption cores 
in spectral lines, as evidence of a grazing eclipse of the accretion disk by the secondary star 
and adopt an orbital inclination of i = 65°. 

Figure 5 shows theoretical radial velocity curves for both cases. The solid curves are the 
radial velocity curves for the q = 0.30 case; the smaller amplitude curve represents the WD 
while the larger amplitude curve is the secondary star. The dashed curve represents the WD 
for the q = 0.20 case, with i = 67°. Gamma velocities of -105 km sec~^ and 8 km sec~^ were 
subtracted from Ha and H/3, respectively. These velocities were determined empirically by 
minimizing the sum of the residuals from the theoretical radial velocity curve. The curves 
show the insensitivity of the system to the value of q. We stress that this paper does not 
argue that the Hct and H/3 radial velocities track the WD motion in SW Sex systems in 
general and in this system in particular. 



Figure 6 shows a projection view of the q 
is closely similar to Figure 6. 



0.30 system. A plot for the q = 0.20 case 
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Panei. Althaos. fc Benvenutd (I2OOOI ) determine a radius of r„d = O.OO771i?0 for a WD 
mass of I.OOM0 and a homogeneous Hamada-Salpeter carbon model. For a secondary star 
mass of O.3OM0, and assuming the secondary is similar to a maii i sequence star of th e same 
mass, we adopt a polar Teg of 3500K for the secondary star (IBaraffe et al.l Il998l ). The 
following analysis will show that the secondary star makes a negligible contribution to the 
system shortward of 7000A. 

Seve ral SW Sex stars have high WD T,^ values. Examples include DW U Ma, 40,000K 



beve rai aw aex stars nave nign wu ipff values, .bxampies mcmae J JW u ivia, 4U,UUUJ\- 
50,000K knigge et allbood krauio-Betancor et al.lboosi : iHoard et allboosh . and UUAqr 



CV, has a Teg of 47,000K flHoard et al 



34,000K ( Baptista. Steiner. fc Cielinski I 1994r) . Similarly, MV Lyr, a well-known novalike 



2004h . We chose a test Tes of 35,000K for SDSSJ0809, 



with the possibility in mind of a subsequent revision upward (or downward). 



Initial estimate of mass transfer rate and the accretion disk tidal cutoff 

radius 



The spectroscopic similarity of SDSSJ0809 to an SW Sex system a l so suggests our initial 

choice of mass transfer rate . SW Sex stars (for reviews see lHellierl (120001 ) : lGroot. Rutten &: van Paradijs 
(I2OOOI ) : iHoard et al.l (l2003l )) are n ovalike (NL) object s with mass tr ansfer rates above the crit- 
ical rate to avoid disk outbursts (Osaki 19961 ). From Osaki ( 1996 ). the critical mass transfer 
rate is Mcrit=4.28x 10-9(Porb/4^'') Meyr^^ In the case of SDSSJ0809, the Osaki expression 
gives a critical rate of 3.4x lO'^ MfD/yr"-*^, but we expect there may be appreciable case-to- 



case variation. 



PattersonI (11984r) predicts a r aass t ransfer rate of 2.5x10 ^^MQyr ^ at the 
top of the period gap, while iHameury et al.l (119881 ) determine a value of l.OxlO~^M0yr"^. 



We discuss annulus models for four mass transfer rates in the following sections. 



A number of studies have considered the tidal cutoff 



j ound a ry, rn, of accretion disks 



( Paczvnski 


(1977 


); 


Papaloizou 


& Pringle ( 


1977 




Whitehurst 


(1988h: 


Schwarzenberg-Czernv & Rozvczka 


(1988 


); 


Whitehurst & King 


( 


1991 


); 


Goodman 




1993 


); 


Warner 


( 


1995)). Equation 2.61 in 



Warned (Il995l ). gives a tidal cutoff radius of 0.49D, where D is the separation of the stellar 
components. The other authors cited agree on rd=0.33D. The latter prescription gives a 
tidal cutoff radius of 51rwd, where r„d is the WD radius. In this connection, the difference 
in D between a g = 0.30 system and a g = 0.20 system is relatively small. Consequently, the 
cutoff radius of the accretion disk is insensitive to q. 

CVs whose mass transfer rates are above the critical rate, and so are s table, are expected 
to ha ve an acc r etion disk temperature profile given by the standard model (jFrank. King, fc Raine 
I992I ). ISmakI (119821 ) finds that the critical temperature for instability in the DIM model 
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( lLasotall200ll ) occurs at Tcfj=6300K. The standard model with the = 0.49D cutoff radius 
produces outer annuh Tea values which are below the critical temperature for instability. 
The = 0.33D cutoff radius produces an outermost annulus T^s value that is at or above 
the critical temperature for a r nass t ransfer rate of at least l.OxlO~^M0yr~^. This result 
agrees with the study by ISmakI (Il983l ). that NL systems have outer radii whose T^s values he 
above the critical temperature. (We thank J. -P. Lasota for a communication on this point.) 
We adopt the tidal cutoff radius = 0.33D for this system, and note that the WD mass 
and system mass ratio are uncertain, with a corresponding uncertainty in the exact cutoff 
radius. 



10. Annulus models for four mass transfer rates 



Without independent information on the mass transfer rate, we used rates typical for 
systems with periods between 3-4 hrs. To this end we have calculated standard models for 
a range of mass transfer rates characteristic of SW Sex stars and have found the best-fitting 
model for each mass transfer rate. We calculated four sets of accretion disk annuli with the 
adopted WD mass of M = l.OM©. We chose mass transfer rates of M=l.OxlO~^M0yr~^, 
2.OxlO-^M0yr-\ 3.OxlO"^M0yr-^ and S.OxlO^^Moyr-^ Table 2 lists annulus properties 
for the mass transfer rate we ultimately chose as most likely. 

The TLUSTY default value for viscosity is a = 0.1, and since this is the value widely 



used for the "hot branch" in the DIM, we adopt it flPringld Il98ll : ISmakI l200d iToutI 12000 



Lasotall200ll ). The viscosity prescripti on strongly affects the annulus verti cal structure, in- 
cluding the annulus optical thickness (jKffz fc Hubenylll986l : iHubenylllQQOl ). The TLUSTY 
parameters Co and (i control the vertical viscosity profile in an a nnulus. We have used the de- 
fault values which assure that there is no "thermal catastrophe" (lHubenylll990l : iHubeny fc Hubeny 
19981 ). The "thermal catastrophe" results when the vertical viscosity profile requires appre- 
ciable energy dissipation in low density layers where cooling from strong resonance lines of 
light metals becomes important. 



The TLUSTY models all produce a standard model Teg (IFrank. King, fc Raindll992l ) 
for an individual annulus. The standard model requires all of the thermal energy gener- 
ated within an annulus to result from viscous dissipation. Our current models also include 
irradiation of the disk by the WD, which can modify the standard model Teff(R) relation. 

The chemical composition of the inner annuli, extending to a transition radius described 
below, include H and He as explicit atoms, and the remaining first 30 atomic species as 
implicit atoms. Inclusion as an implicit atom means that the atomic species contributes 
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to the total number of particles and the total charge, but not to the opacity. The vertical 
temperature profile of a given small radius annulus varies smoothly from the central plane to 
the upper boundary of the annulus. At larger radii the temperature profiles show increasingly 
steep drops toward the upper boundary. For these annuli we "switched on" convection in 
TLUSTY. In addition, the larger radius annuli have lower T^s values, and opacity due to 
metals is included. Beginning at Teff~12000K, our annulus synthetic spectra include H, He, 
C, Mg, Al, Si, and Fe as exphcit atoms and the remaining first 30 atomic species as implicit 
atoms. Because they have fairly high-lying first excited states, N and O do not contribute 
importantly to the continuum opacity of interest (photoionization from the ground state is 
close to the Lyman limit, indeed shortward of it for O). 

The data for a given annulus have been extracted from an extensive TLUSTY tabulation 
of physical properties at all (typically 70) reference levels in the annulus, from the central 
plane to the upper boundary. Note from the last column of Table 2 that all annuli are 
optically thick. This condition, that all annuli are optically thick, is also true for all of the 
other mass transfer rates calculated. 

The outer annuli for the mass transfer rate of LOxlO^^Mgyr"^ have outermost annulus 
Teff values that fall below the stabihty hmit for the DIM model. Smaller mass transfer rates 
have appreciably larger fractions of the outer annuli that fall below the stability limit. This 
sets a lower limit on the mass transfer rates to consider. 



11. System models for successive mass transfer rates 

We began with a standard model T{R) for a mass transfer rate of LOxlO^^M© yr~^ 
and an untruncated accretion disk. Because the FUSE2 exposure was closer to the time 
of the HST exposure than the FUSEl exposure, we initially chose the FUSE2 spectrum 
for comparison with our simulations. (The number of plots required to show the results for 
all four mass transfer rates would be very large. Consequently we present illustrative plots 
for only one of the mass transfer rates. Similarly, we tabulate annulus properties for only 
one of the four mass transfer rates (Table 2)). A continuum spectrum fit to the FUSE2, 
HST and SDSS spectra for the untruncated disk had far too much UV flux as compared 
with the HST and FUSE2 spectra. This implied too large a contribution from high T^s 
annuli. From equation (1), the only way to reduce those contributions, and maintain the 
standard model, is to truncate the inner accretion disk. Accordingly, we started with the 
outer radius at the tidal cutoff radius and calculated models truncated at inner radii of 1.0, 
4.0, 6.0, 7.5, 8.2, 9.0, and 10.0 times the WD radius. We calculated continuum spectra for the 
models and determined individual scaling factors to fit them to the SDSS spectrum, requiring 
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accurate superposition among the models at the long wavelength end of the SDSS spectrum. 

A model truncated at 8.2 times the WD radius was the best overall fit to the FUSE2, 
HST, and SDSS spectra, as judged visually. It would be very difficult to be quantitative in 
this estimate, since the degree of departure from a good fit in the optical must be weighed 
against what is an approximate fit to the HST and FU SE spectra. In addition, the many 
UV emission line features make location of the continuum problematic. The illustrative plots 
explain this difficulty. The synthetic spectrum spectral gradient was clearly too shallow in 
the fit to the SDSS spectrum, and was shghtly too hot in the UV. The UV discrepancy 
increased on substituting the FUSEl spectrum. Since the FUSEl spectrum provides a 
more stringent UV constraint on the models (setting a maximum WD contribution), and 
since the FUSEl spectrum mates fairly smoothly with the HST spectrum in their overlap 
region, while the FUSE2 spectrum docs not, we decided to use the FUSEl spectrum in 
our final comparisons between the models and the observed spectra. It is possible that 
the higher excitation FUSE2 spectrum results from a higher mass transfer rate from the 
secondary star. 

A slight increase in the truncation radius, to 9.0 times the WD radius, made the fit to 
the optical spectrum slightly worse without a great improvement in the UV fit. Since the 
system synthetic spectrum had too much UV flux, and since increasing the truncation radius 
did not help, we tested the other option of reducing the WD T^g- Substituting a 20,000K 
WD reduced the WD contribution to the synthetic spectrum to a nearly negligible value, 
and improved the UV fit. 

A mass transfer rate of 2.0 x lO^^M© yr^^ presented the same problem for an untruncated 
accretion disk. Following the procedure of the previous mass transfer rate, we started with 
the outer radius at the tidal cutoff radius and calculated models truncated at inner radii of 
1.0, 4.0, 6.0, 8.0, 9.7, and 12.0 times the WD radius. The 9.7 model produced a very good fit 
to the optical (SDSS) spectrum but was too hot in the UV. As with the first mass transfer 
rate, substituting a 20,000K WD substantially improved the UV fit. 

We next considered a mass transfer rate of S.OxIO^^Mq yr~^. We started with the 
outer radius at the tidal cutoff radius and calculated models truncated at inner radii of 1.0, 
4.0, 8.0, 12.0, 14.3, 16.0, and 18.0 times the WD radius. Figiuc 7 shows the superposition 
of these models on the optical spectrum. Note the progression from a too large spectral 
gradient to a too small gradient for the synthetic spectra. It clearly is possible to select 
a best fitting model in this spectral region, by visual inspection, from among the available 
choices. Although we have a fitting procedure available, we do not believe its use is 
warranted in the present illustration. Figure 8 connects the optical region to the UV. Note 
the large separation of the synthetic spectra in the UV for relatively small differences in the 
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optical. Figure 9 shows the fits in the UV only. The best overall fit appears to be the model 
with inner truncation radius r = IG.Or^d- (For brevity, the subsequent text will use "the 
16.0 model" to reference this model, and similarly for other inner truncation radii.) The 
14.3 model is too bright in the llOOA to 1350A interval while the 18.0 model has a large 
flux discrepancy longward of 1250A. Figure 10 shows the fit of the 16.0 model to the optical 
spectrum, and Figure 11 shows the fit to the UV spectra. The Figure 10 result, that the 
contribution of the secondary star is essentially negligible in the optical, and completely so 
in the UV, applies to all of the models in this investigation. The fit of the 14.3 model to the 
optical spectrum is nearly identical to Figure 10. From Figure 7, the 18.0 model (the lowest 
of the seven plots) has a clearly too small spectral gradient in the optical. This discussion 
illustrates how we can use visual estimates to select the best fitting model. Plots for other 
mass transfer rates show characteristics comparable to those shown here but with specific 
differences described for those mass transfer rates. 

For comparison. Figure 12 substitutes the FUSE2 spectrum. Note that, in Figure 11, 
the FU SEl spectrum and the HST spectrum fit reasonably smoothly in the overlap region, 
1150A to 1200A, while there is a clear discontinuity in the Figure 12 overlap. This comparison 
justifies our use of the FUSEl spectrum. An uncertainty remains for the SDSS spectrum 
because of its non-simultaneity with the HST or FU SE spectra. 

The fit to the FUSEl spectrum shown in Figure 11 places the synthetic spectrum too 
high in the vicinity of 1150A. Figure 13 shows the result of substituting a 20,000K WD. 
The WD contribution, barely visible at the bottom, now is nearly negligible. The fit near 
1150A is improved but the fiux deficit longward of 1300A is larger. 

The same procedure was followed for the mass transfer rate of 5. Ox 10^*^ yr^^. We 
calculated a set of five models which bracketed the HST spectrum from too much to too 
little fiux, with the synthetic spectra normalized to agree with the optical spectrum at 
9200A. All of these models showed a too large spectral gradient on comparison with the 
optical spectrum. All of the models were too hot in the FUSE range, for a 35,000K WD. 
Substituting a 20,000K WD improved the UV fit but in the FUSE range the synthetic 
spectrum still was too hot. With this substitution, as with the other mass transfer rates, 
the WD contribution became neghgible. 

Larger mass transfer rates would produce still greater discrepancies in the fits. Con- 
sequently, we do not consider larger mass transfer rates. The best-fitting model at the 
minimum mass transfer rate, l.OxlO~^M0 yr~^, is clearly inferior to the best-fitting models 
at 2.OxlO~^M0 yr~^ and S.OxlO^^M© yr^^, which are nearly comparable. The best-fitting 
model at the maximum mass transfer rate, S.OxlO^^M© yr~^, also is inferior to the two 
intermediate mass transfer rates. We judge that the 3.OxlO~^M0 yr~^ mass transfer rate 
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is a slightly better spectral fit than the 2.0x10 ^Mq yr ^ rate and adopt the mass transfer 
rate of S.OxIO^^Mq yr"-*^ for this system. 

The most important result from the analysis of the four mass transfer rates is that 
the accretion disk must be truncated if the standard model applies. The fits to the HST 
spectrum exhibit a flux deficiency between about 1400A and the 1700A limit of the HST 
spectrum for models that have acceptable flux levels in the FUSE range. We see no way 
to accommodate that discrepancy by maintaining the standard model. If this result is to 



possibilities are evaporation (IMenou et al.l (120001 ): iDubus et al.l (I2OOII ): 



Lasota 


(2001 


)) and 


Menou et al.l 


(2000. 



interaction with a WD magnetic field. Using the evaporation equation of 
eq.3), and setting the evaporation rate equal to the mass tra nsfer rate, it can be shown that 



evapo ration is too small an effect by four orders of magnitude. iGroot. Rutten fc van Paradijs 



( I2OOOI ) sta te that the SW Sex phenorn enon can be explained in the context of a non-magnetic 



WD, and iRodriguez-Gil et al.l (120011 ) identify circular polarization in the SW Sex star LS 



Peg, and th e y arg ue that magnetic accretion plays a fundamental role in SW Sex stars. 



Hoard et al.l (l2003l ) summarize and describe the various arguments and evidence for and 
against magnetic WDs in SW Sex stars. Because of this uncertainty, and in anticipation of 
the results presented in §12, we do not discuss possible truncation of the accretion disk by a 
magnetic field. 

A further important result from the truncation models is that the best fit to the FU SE 
spectrum requires a nearly negligible WD contribution to the system synthetic spectrum. 
This result is hard to understand, since Figure 6 clearly shows that the inclination is too 
small for a vertically enlarged accretion disk rim to hide the WD. Hiding the WD would 
be even more difficult if the accretion disk is truncated. The negligible WD contribution 
cannot be a result of the accretion disk hiding the WD. The truncation models cover the 
full range of acceptable mass transfer rates, and a cool WD produces a better synthetic 
system spectrum fit than the 35,000K WD in all cases. The alternative to hiding the WD 
behind the accretion disk rim, that the WD actually is cool, is difficult to under st ar id in 
view of the large raass t r ansfer rate which would be e xpected to produce a hot WD (ISion 



I999I: ISzkodv et al 



2002 



Arauio-Betancor et al.ll2003l ). The spreading layer theory applied 



to WDs (jPiro &: BildstenI 12004 ) might permit a prediction of the WD T^g for given Mwd and 
M. Although we have found a model which is a fairly good fit to the combined spectroscopic 
data, the model appears unrealistic and we reject the truncated accretion disk scenario. 

We are forced to conclude that the accretion disk in SDSSJ0809 may depart from a 
standard model, which is the basis of all the models we have calculated to this point. 



- 18 - 



12. Non-standard model accretion disks 



Various studies (e.g., Rutten et al. ( 19921 )). based on image reconstruction, find a flatter 



temperature proflle than the standard T{R) relation. It is of interest to adopt an isothermal 
accretion disk and test whether that prescription produces a viable model. Figure 14 shows 
the flts of isothermal 15,500K, 13,000K, and 12,000K untruncated accretion disks to the 
observed optical spectrum. Figure 15 shows the flts to the HST and FUSEl spectra. Each 
model includes the contribution of a 35,000K WD. The WD contribution, normalized as 
for the 13000K model, is shown separately at the bottom of Figure 14 and Figure 15. The 
12,000K model flts the HST spectrum fairly well but has too little flux in the FU SE region. 
All three models are poor flts to the observed optical spectrum. 



The tomographic image reconstruction for SW Sex itself (iRutten et al.l (Il992l )) shows 
that the accretion disk temperature proflle follows the standard model in the outer part 
of the accretion disk, but at a transition radius the temperature proflle becomes flat and 
remains nearly so to the inner edge, with a slow increase. We have implemented a model 
similar to this reconstruction by flrst adopting the mass transfer rate, 3.0xl0~^ Mq yr~^, for 
the best standard model, and then flnding the radius for which the local T^s was 13,000K. 
This transition value was based on our result for the isothermal disk models. We used an 
accretion disk that was not truncated and set the T^s of all annuli within the transition radius 
equal to 13,000K. The resulting model, in contrast to the isothermal models, flts the optical 
spectrum very well. However, the flt to the HST spectrum left appreciable positive residuals. 
We experimented with a larger temperature isothermal region, with a smaller radius for the 
crossover to the standard model. Eventually we found that an isothermal region of 14,000K, 
with a crossover to the standard model at r/r„d = 18.0, gave a fair representation of the 
HST and FUSEl spectra. (Note that the tidal cutoff radius is at r/r^d = 51.3). The 
temperature proflle of the accretion disk is in Table 5. This model assumes the accretion 
disk extends to the WD equator, so only half of the WD is visible. The BINSYN software 
allows for eclipse of inner annuli segments that the WD eclipses. 

The flt to the HST and FUSEl spectra suggest the possibility of an improved flt with 
a higher T^s WD. We calculated models, in addition to the 35,000K model, of 40,000K, 
45,000K, and 50,000K. All four models have nearly identical flts to the upper (discovery) 
SDSS spectrum, as shown in Figure 16. The lower (archival) SDSS spectrum is discussed 
in the following section. Figure 17 shows the flts to the HST and FUSEl spectra. The 
bottom spectrum is the contribution of the accretion disk. Based on the flt near 900A, the 
45,000K model roughly bisects the observed continuum. 

The synthetic spectra are given in Eddington flux units; a divisor, sn, is applied to 
superpose a synthetic spectrum on the observed spectra. Since the observed flux is tabulated 
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in units of 1.0x10 erg s ^ cm ^ A ^, the distance, d, is given by (P — sat/I-OxIO 

In the present case, for the model using the discovery SDSS spectrum, the scahng factor is 
6.1x10^8; the formal derived distance to SDSSJ0809 is 800pc. Note that there is a large 
associated uncertainty. Tentative system parameters are listed in Table 3. The Table 3 
parameters are consistent with other SW Sex systems. 

13. The archival SDSS spectrum and model implications 

Figure 16 shows the fit of the Table 5 model to the discovery SDSS spectrum. The 
lower SDSS spectrum is the archival spectrum. Dividing the discovery spectrum by 1.3 
and overplotting produces a close fit to the archival spectrum between 5500A and 9200A. 
The discovery spectrum shows increasing positive residuals shortward of 5500A, reaching a 
maximum at 4100A, and decreasing to coincidence with the archival spectrum at 3800A. 

Before attempting to model the archival SDSS spectrum we seek a likely physical basis 
for the observed change from the discovery spectrum. An important point is that the change 
took place in only two days. In §11 we have argued that the accretion disk cannot be approx- 
imated by a standard model and that models truncated on the inner edge are unrealistic. In 
§12 we found that isothermal accretion disks produce poorly-fitting models, but that a model 
whose temperature profile mimics SW Sex satisfactorily fits the discovery SDSS spectrum 
as well as the UV data. In §11 we have already suggested that the FUSE2 spectrum differs 
from the FUSEl spectrum because of a higher mass transer rate. We believe the two SDSS 
spectra most likely differ for the same reason. 

If the mass transfer rate from the secondary star decreases, a cooling wave will propagate 
inward from the outer accretion disk rim. The temperature of the outer annuli could fall 
below the transition temperature from the high state. With this scenario in mind, we started 
with the Table 5 profile, reduced the outer annuli temperatures and calculated a new system 
model. Although the fiux level was reduced, the system synthetic spectrum had a spectral 
gradient that was too high compared with the archival spectrum. The scaling factor to 
superpose the synthetic spectrum on either observed spectrum has to remain fixed, since 
it connects directly to the system distance. We continued temperature reductions for all 
annuh, reaching the final temperature profile of Table 6. Note that a complete system model 
was required at each new accretion disk temperature profile. The final fit to the archival 
spectrum is in Figure 18. Note that Figure 18 presents fits for the same four values of the 
WD temperature as in Figure 16. Figure 19 shows the four model fits to the UV spectra. 
Note the large reduction in the contribution of the accretion disk as shown by the lowest 
plot. 
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These fits have unacceptably large residuals from the HST spectrum. A basic problem 
is that we have only a single HST spectrum and it is not contemporaneous with any of 
the other observed spectra. No information is available about whether the HST spectrum 
changes in a time interval of days, and by how much. 

The proposed scenario to explain the differences in the two SPSS spectra resembles the 



widely-accepted explanation of VY Sculptoris star behavior (ILivio fc Pringlelll994j : iKing fc Cannizzo 



19981 ). However, SW Sex objects do not undergo the very large brightness variations of 
VY Scl systems, so the physical cause of mass transfer variation in SDSSJ0809 is un- 
certain. A detailed study of the high-state/ l ow-st ate transition in the VY Scl system 



V794 Aql (IHoneycutt. Cannizzo. &: Robertson! Il994l ) provides a useful comparison. The 



drops in V794 Aql brightness have rates that exceed the two-day drop of 0.3 mag. for 
SDSSJ0809, so the SDSSJ0809 drop rate is not excessive. The theoretical time-dependent 
model for V794 Aql shows two-day changes in annulus T^s values as large as 10,000K. Also, 
successive two-day log T^s profiles across the accretion disk are nearly vertical displacements 
from each other except near the end of the simulation. Thus, the changes from Table 5 to 
Table 6 appear to be theoretically acceptable, based on a comparison with V794 Aql. The 
two cases are not entirely comparable since the SDSSJ0809 starting model (Table 5) is not 
a standard model. 

Although it originally appeared reasonable to calculate a system model fitting the com- 
bined spectra then available, we conclude that contemporaneous spectra are essential for 
reliable models of systems which may show significant day-to-day changes. Note that this 
conclusion is based on a quantitative comparison with model calculations. 



14. Discussion 

The FUSEl spectrum fits smoothly with the HST spectrum in the overlap region, 
while the FUSE2 spectrum does not fit well. This is inconclusive evidence that SDSSJ0809 
was in a lower accretion state at the times of both the HST and FUSEl exposures, and that 
the FU SE2 spectrum was obtained during a higher state. Because of this uncertainty, we 
emphasize again that this study has been an illustration of techniques to simulate accretion 
disks, with the SDSSJ0809 observational data serving as the vehicle for that illustration. 

All of these TLUSTY annulus models are optically thick, including the nonstandard 
models. The optically-thick annulus spectra with their absorption lines and large Balmer 
discontinuities are in conflict with all of the observed emission line spectra {FUSE, STIS, 
and SDSS). What do the emission lines tell us and what is the most likely way the model 
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should be changed to produce the observed emission hues? All of the SDSSJ0809 emission 
lines are single-peaked, one of the defining characteristics of SW Sex stars. In contrast to 
the SW Sex result that the emission l ines arise excl usively from the accretion stream shock 
( IGroot. Rutten fc van Paradijsll2000l ). iHellierl (120001 ) argues for emission lines produced in a 
wind. Hellier asserts that a wind provides the most straightforward explanation of single- 
peaked emission lines, and that in particular the wind component fills in any double-peaked 
disk emission in the Balmer lines. (Also see iDrew fc Progal (120001 )). Table 4 provides some 
details concerning the Balmer emission lines in SDSSJ0809. There are no P Cygni profiles 
among the SDSSJ0809 Balmer lines or He lines. There is only a hint of an absorption 
component at the blue edge of the C IV 1550A lines. 

The most direct explanation of the SDSSJ0809 emission lines is production in a chro- 
mospheric or coronal region overlyi ng the accretion disk. Severa l mech anisms have been pro- 
posed that cou l d pro duce coronae: iMeyer fc Meyer-Hofmeisteii (119941 ) use a coronal siphon; 
Murray fc Lin] (119921 ) depend on sound waves that accelerate to form shocks; iLiang fc Price 
( 1l977l ) pr opose general n onthermal processes to transport energy vertically from the accre- 
tion disk: iHubenvi (1l990l ) depends on the vertical variatio n of viscosity within an annulus; 
Cannizzol J2000h applies a thermal evaporative instability (jShaviv fc Wehrsdll986l ). 



An extensi ve literatu r e exis ts rel ating to forrn a tion o f emission lines in CV systems. 
Wilhamsl Jl98oh : iTvlendal Jl98lh and Icheng fc LiJ Jl989l) discuss emission line formation 
in optically thin outer annuli of accretion disks. iKo et al.l (119961 ) discuss emission line for- 



mation as a resu 



t of X-ra y /EUV illuminat ion of the accretion disk by the central object. 



Horne &: Marsh] (119861 ) and iFergusonI (119971 ) discuss emission line formation in more general 



cases. The SDSS spectrum does not extend far enough in the UV to determine the last 
resolvable Balmer line; hence, it is not possible to apply the Inglis- Teller equation to deter- 
mine the electron density where the continuum is produced. Note from Table 4 that the 
Balmer decreme nt is relatively fiat. (Compare with case A and case B for planetary neb- 
ulae. Table 2 of lAUer &: Lillerl (jl968l ). and the result that Balmer de c remen ts for some Be 
stars appr oximate those of planet ary nebulae (IBurbidge fc Burbidgd (119531 ) and references 
therein).) iHorne fc MarshI (119861 ) point out that in the case of a fiat Balmer decrement 
the Balmer emission lines are optically thick. The fact that all of the annulus models are 
optically thick, together with the result that the Balmer emission lines are optically thick, 
is consistent with an origin of emission lines in a chromosphere or corona, not in optically 
thin outer annuli. Moreover, emission lines produced in optically thi n outer regions of a n 
accretion disk occur in a temperature environment of 6000K-9000K (jCheng fc Linlll989l ). 
The high excitation emission lines in SDSSJ0809, particularly in the FUSE spectra, require 
a much higher temperature. A notable feature of the observed optical Balmer emission lines 
is that they are narrower than the synthetic spectrum absorption lines (e.g., see Figure 10). 
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The model absorption lines are broadened by Keplerian rotation appropriate to the orbital 
plane. The narrower emission lines could be an indication they are produced well above the 
accretion disk face, where the horizontal acceleration due to the WD, and the associated 
Keplerian speed, is smaller. 



Extrapolation of the SDSS spectrum indicates the absence of a Balmer jump. Ila Dous 



(Il989l ) points out that this condition implies a much lower vertical temperature gradient 
in an accretion disk than is true in a stellar atmosphere, and likely associates with energy 
production in the accretion disk, in contrast to a stellar atmosphere. 

We doubt that addition of an emission line slab to any of the models will smoothly fill 
in the absorption lines and convert the absorption line spectrum to the observed emission 
line spectrum. Rather, we believe that a modification of the individual annulus models is 
neces sary to convert their spectra to emission line spectra. What are the prospects to do 
this? iHubenyl ( 1l990l ) provides an analytical model for accretion disk annuli that is the basis 
of the TLUSTY annulus model calculation. As discussed in §10, TLUSTY includes two 
parameters, Co and Ci, that specify the vertical viscosity profile in an annulus, and therefore 
the rate of energy generation as a funct ion of z value. These parameters have de fault values 
that prevent a "thermal catastrophe" (iHubenyl Il990l : iHubeny fc Hubenyl 119981 ). Carefully 
modified values of these parameters can, in principle, produce emission lines. 



The IHubenyl ( Il990l ) analytical model specifically considers co nditions t hat le ad to bound- 
ary temperatures higher than the annulus Tgfj. As discussed by IHubenyl ( Il990l ). the energy 
dissipated in low density (coronal) layers can be removed by strong (emission) resonance 
lines of abundant species, such as H I, Mg II, C II, Al HI, Si IV, C IV, N V, O VI, etc. 
Numerical simulations for actual physical models (IStone et al.lll996l ) indicate t hat vi s cosity 
increases toward the surface, as required for the prop osed TLUSTY models. iToutI (120001 ) 
discusses our present understanding of viscosity, while iHawleyl (120011 ) discusses the present 
status of magnetohydrodynamic simulations of cylindrical Keplerian disks and their relation 
to viscosity. The calculation of TLUSTY models with modified values of Co and Ci is beyond 
the scope of the present investigation. 

We have not accounted for the half of the liberated potential energy associated with 
the mass transfer stream that is usually assigned to accretion luminosity. Part could be 
dissipated in the shock where the transfer stream merges with the accretion disk. Although 
BINSYN has the capability to model a rim bright spot, the light curve provides no indication 
of a bright spot's existence. A major part of the energy goes into heating the WD. We assume 
the remainder is carried away in a wind. 



An original motivation for analyzing SDSSJ0809 was the hope that visibility of the 
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WD could be demonstrated. The primary argument for detection of the WD is the model's 
approximate representation of the UV continuum shortward of 1200A. If more extensive 
observational data support this model, then SDSSJ0809 is very likely an SW Sex star with 
a high WD Teg in agreement with T^g values for other SW Sex stars. 
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grateful for support from NASA FUSE giant NAG5-13656, NASA HST grants GO-09357.06, 
GO-09724, and AR-10674, and NSF grant AST 02-05875. BTG was supported by a PPARC 
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Table 1. Observation Log 



Observation 


Exposure 


Start Time 


Total Exposure 


# 


# 


(UT) 


(HJD-2450000) 


(s) 


1 


1 


28 Mar 2003 11:36:10 


2726.985218 


1590 




2 


28 Mar 2003 13:01:23 


2727.044391 


2473 




3 


28 Mar 2003 14:40:37 


2727.113297 


2515 




4 


28 Mar 2003 16:29:19 


2727.188776 


1888 


2 


1 


16 Mar 2004 14:44:13 


3081.116762 


3257 




2 


16 Mar 2004 16:44:41 


3081.200413 


2023 




3 


16 Mar 2004 18:31:30 


3081.274585 


1606 




4 


16 Mar 2004 20:09:14 


3081.342450 


1737 




5 


16 Mar 2004 21:53:07 


3081.414585 


1495 




6 


16 Mar 2004 23:02:56 


3081.463065 


2050 




7 


16 Mar 2004 23:36:45 


3081.486547 


1272 




8 


17 Mar 2004 00:43:58 


3081.533221 


560 



Note. — "Total Exposure" will be longer than the final usable exposure time because 
of data that are rejected for quality reasons during the standard pipeline processing. 
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Table 2. Properties of accretion disk with mass transfer rate M — 3.0£'(— 9)MQyr ^ and 

WD mass of l.OMp,. 









mo 


To 


log 


g 




zo 




Ne 


''"Ross 


1, 


.36 


66333 


8. 


,967E3 


54213 


7, 


.16 


4. 


,63E7 


1, 


.11E15 


1. 


,48E4 


2, 


.00 


59475 


1. 


,095E4 


48448 


6, 


.90 


8. 


,01E7 


6, 


.75E14 


1. 


,76E4 


3, 


.00 


48093 


1. 


,070E4 


39159 


6, 


.60 


1. 


,35E8 


4, 


.14E14 


2. 


,08E4 


4, 


.00 


40423 


9. 


,959E3 


32916 


6, 


.38 


1. 


,94E8 


2, 


.90E14 


2. 


,39E4 


5, 


.00 


35062 


9. 


.251E3 


28556 


6, 


.21 


2. 


.55E8 


2, 


.24E14 


2. 


.71E4 


6, 


.00 


31106 


8. 


.639E3 


25342 


6, 


.07 


3. 


.18E8 


1, 


.82E14 


3. 


.06E4 


7, 


.00 


28058 


8. 


.115E3 


22867 


5, 


.95 


3. 


.81E8 


1, 


.52E14 


3. 


.39E4 


8, 


.00 


25630 


7. 


.665E3 


20894 


5, 


.84 


4. 


.47E8 


1 


.31E14 


3. 


.77E4 


1 n 

lU, 


nn 
.uu 


21986 


6. 


,934E3 


17395 


5, 


.67 


5. 


,84E8 


1, 


.02E14 


A 

4. 




12, 


.00 


19367 


6. 


,365E3 


15812 


5, 


.53 


7, 


,24E8 


8, 


.21E13 


5. 


,42E4 


14, 


.00 


17381 


5. 


,907E3 


14201 


5, 


.40 


8, 


,68E8 


6, 


.74E13 


6. 


,25E4 


16, 


.00 


15816 


5, 


,529E3 


12933 


5, 


.30 


1. 


,01E9 


5, 


.71E13 


7, 


,00E4 


18, 


.00 


14547 


5. 


,210E3 


11905 


5, 


.20 


1. 


,16E9 


4, 


.96E13 


7. 


,60E4 


20, 


.00 


13495 


4. 


.937E3 


11054 


5, 


.12 


1. 


.31E9 


4, 


.36E13 


7. 


.98E4 


22, 


.00 


12606 


4. 


.700E3 


10333 


5, 


.04 


1. 


.44E9 


3 


.78E13 


8. 


.13E4 


24, 


.00 


11843 


4. 


.492E3 


9717 


4, 


.97 


1. 


.59E9 


3 


.25E13 


8. 


.10E4 


30, 


.00 


10086 


3. 


,993E3 


8287 


4, 


.78 


2. 


,02E9 


1, 


.49E13 


7. 


,52E4 


34, 


.00 


9213 


3. 


,734E3 


7572 


4, 


.68 


2, 


,32E9 


6, 


.68E12 


7. 


,17E4 


40, 


.00 


8188 


3, 


,420E3 


6744 


4, 


.61 


3. 


,13E9 


1, 


.87E12 


6, 


,97E4 


45, 


.00 


7516 


3, 


,207E3 


6202 


4, 


.39 


2. 


,68E9 


5, 


.27E11 


6. 


,69E4 


50, 


.00 


6961 


3. 


,026E3 


5759 


4, 


.18 


2. 


,36E9 


1, 


.59E11 


7. 


,48E4 



Note. — Each line in the table represents a separate annulus. The 
column headed by mo is the column mass above the central plane. The 
column headed by To is the boundary temperature. Compare with Teff- 
The log g values are at an optical depth of 0.9. The zq column gives the 
height of the anmihis in cm. The accretion disk radius at the tidal cutoff 
boundary is 2.75xl0^°cm. The Ne column is the electron density at the 
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upper annulus boundary. The tross column is the Rosseland optical 
depth at the central plane. A viscosity parameter a = 0.1 was used for 
all annuli. 
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Table 3. SDSSJ0809 Model System Parameters 



parameter value parameter value 





l.OO±O.2M0 


''"wd 


0.00771i?( 




0.30±0.10Mq 


log fi-wd 


8.35 


M 


3.O±l.OxlO-''M0yr-i 


rs(pole) 


0.3117^0 


P 


0.133 days 


rs(point) 


O.451i?0 


D 


1.19Oli?0 


rs(side) 


O.324i?0 




155.7 


Tsback) 


O.363i?0 




2.46623 


log (pole) 


4.95 


i 


65±5° 


log (point) 


-0.44 




45,000±5000K 


log 5's(side) 


4.88 


Teff,s(pole) 


3500K (nominal) 


log ^^(back) 


4.67 


^wd 


1.0 




O.39i?0 


As 


0.5 


n 


0.00771i?( 


bwd 


0.25 


H 


O.OO96i?0 


hs 


0.08 







Note. — wd refers to the WD; s refers to the secondary star. 
D is the component separation of centers, f2 is a Roche potential. 
Temperatures are polar values, A values are bolometric albedos, 
and b values are gravity-darkening exponents, specifics the outer 
radius of the accretion disk, set at the tidal cut-off radius, and Vh is 
the accretion disk inner radius, as determined in the final system 
model, while H is the semi-height of the accretion disk rim, based 
on the standard model. 



-32- 



Table 4. Features of SDSSJ0809 Balmer emission lines 



Identification 


FWHM 


EW 


decrement 


Ha 


21.5 


21.6 


2.5 


H/3 


14.0 


8.78 


1.00 


H-f 


13.0 


5.98 


0.68 


H5 


12.0 


4.85 


0.55 


He 


14.0^ 






HS 


ll.O'' 






H9 


6.0 


1.00 


0.11 



^Blended with He II 3968.43A 

''Blended with He 11 3887.44A 



Table 5. Accretion disk temperature profile for discovery SDSS spectrum 















1.0000 


14,000 


16.3448 


14,000 


34.6581 


9079 


1.0225 


14,000 


18.0096 


14,000 


36.3230 


8776 


1.3611 


14,000 


19.6745 


13,644 


37.9878 


8495 


3.0260 


14,000 


21.3393 


12,875 


39.6527 


8234 


4.6908 


14,000 


23.0042 


12,200 


41.3175 


7992 


6.3557 


14,000 


24.6690 


11,164 


42.9824 


7766 


8.0205 


14,000 


26.3339 


11,071 


44.6472 


7554 


9.6854 


14,000 


27.9987 


10,593 


46.3121 


7355 


11.3502 


14,000 


29.6636 


10,162 


47.9769 


7168 


13.0151 


14,000 


31.3284 


9768 


49.6418 


6992 


14.6799 


14,000 


32.9933 


9409 


51.3066 


6826 
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Table 6. Accretion disk temperature profile for archival SDSS spectrum 







r/r^d 


TeS 


r/r^d 


TeS 


1 


.0000 


11,500 


16. 


,3448 


10,000 


34. 


,6581 


9000 


1 


.0225 


11,500 


18. 


,0096 


10,000 


36. 


,3230 


8500 


1 


.3611 


11,500 


19. 


,6745 


10,000 


37. 


,9878 


8300 


3, 


.0260 


11,500 


21. 


,3393 


10,000 


39. 


,6527 


8100 


4, 


.6908 


ILOOO 


23. 


,0042 


10,000 


41. 


,3175 


7800 


6, 


.3557 


ILOOO 


24. 


,6690 


10,000 


42. 


,9824 


7500 


8, 


.0205 


11,000 


26. 


,3339 


10,000 


44. 


,6472 


7200 


9, 


.6854 


11,000 


27. 


,9987 


10,000 


46. 


,3121 


7000 


11 


.3502 


10,500 


29. 


,6636 


10,000 


47, 


.9769 


6700 


13 


.0151 


10,500 


31. 


,3284 


9500 


49. 


,6418 


6500 


14, 


.6799 


10,500 


32. 


,9933 


9200 


51. 


,3066 


6200 



(N 

I 



o 

7 

i-H 




900 920 
HI ■ 

Hell 



940 



960 



980 



1000 



J L 



\ L 




1100 



1120 1140 

Wavelength (A) 



1160 



1180 



Fig. 1. — Average FUV spectrum of SDSSJ0809 from Observation 1 (2003), binned to a 
dispersion of 0.20 A pixel"^. The spectrum spans 915-1185 A; the top and bottom panels 
overlap the ends of the middle panel by 7.5 A. Prominent features often found in FUV spectra 
of CVs are indicated, alt hough not all of thern are present in SDSSJ0809 (see text). Airglow 
lines (®; identified from iFeldman et al.ll200ll ) and probable ISM lines (ions in italics) are 
indicated below the spectrum. The airglow lines have been truncated at the upper flux limit 
of the plot. Wavelengths of H I and He II transitions are shown between the panels. 
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Fig. 2. — As in Figured] but showing the average FUV spectrum of SDSSJ0809 from Obser- 
vation 2 (2004). 
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Fig. 3. — i/ST spectrum of SDSSJ0809. Note that several emission features iiave absorption 
cores extending below the continuum, indicating a substantial line-of-sight column density. 
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Fig. 4. — Discovery SDSS spectrum of SDSSJ0809. Note tlie lack of absorption wings on 
the Balmer lines and the absence of any indication of a Balmer discontinuity. 
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Fig. 5. — Theoretical radial velocity curves of system stars compared with Ha (diamonds) 
and up (x's) observations. The WD mass is I.OMq, the mass ratio is q = 0.30, the orbital 
period is 192™, and the orbital inchnation is 65° (continuous curves). The dashed curve is 
the projected WD velocity for a mass ratio of g = 0.20. 
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Fig. 6. — View of the system described in Figure 5 at orbital phase 0.0, q — 0.30, projected 
on the plane of the sky for i — 65°. This is inferior conjunction for the secondary star; it is 
between the observer and the WD. The inner edge of the accretion disk extends to the WD 
equator in this plot (no truncation), and the outer accretion disk radius is at the tidal cutoff 
radius. The plot for a g = 0.20, i = 67° system is closely similar. 
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4000 5000 6000 7000 8000 9000 

Wavelength 

Fig. 7. — Synthetic spectra for truncation radii of (from top to bottom) 1.0, 4.0, 8.0, 12.0, 
14.3, 16.0, and 18.0 times the WD radius (thin lines), fitted to the SDSS spectrum (broad 
fine). The mass transfer rate is 3.0x10~^Mq yr~^. The best fit is the r = 16.0rwD spectrum; 
it shows separately in Figure 10. 
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1000 2000 3000 4000 5000 

Wavelength 

Fig. 8. — As in Figure 7, but extending into the UV region. Note the large separation of 
synthetic spectra in the UV for a relatively small difference in the optical. The mass transfer 
rate is 3.0x10'^ Mq yr'^. 
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Fig. 9. — As in Figure 8 but showing only the UV. The 35,000K WD model spectrum is at 
the bottom. The other model spectra are labeled with the disk inner truncation radius (in 
WD radii). 
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Fig. 10. — Plot showing detail of the Tinner = IG.OrwD model(dark line) fit to the opti- 
cal spectrum (light line). The spectrum rising to the bottom right is the secondary star 
component, and the spectrum rising to the bottom left is the 35,000K WD component. 
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Fig. 11. — As in Figure 10 but for tiie UV data. Note that tiie FUSEl and HST spectra 
agree well in the overlap region. The 35,000K WD contribution is the lowest spectrum. 
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Fig. 12. — As in Figure 11 except the FUSE2 spectrum has been substituted for the 
FUSEl spectrum. The HST and FUSE2 spectra do not match weU in the overlap region. 
The synthetic spectrum is a poor fit to the FU SE2 spectrum. 
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Fig. 13. — As in Figure 11 except a 20,000K WD has been substituted. The synthetic 
spectrum now passes through the middle of the FUSEl spectrum near 1150A, but the flux 
deficit longward of 1300A is appreciable. 



-47- 




4000 5000 6000 7000 8000 9000 

Wavelength 



Fig. 14. — Synthetic spectrum for tliree isotliermal accretion disks (untruncated) . Note that 
the spectral gradient is too large in all cases. 
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Fig. 15. — As in Figure 9 but showing the isothermal disk models from Figure 14. The 
relative contribution of the 35,000K WD is less than in Figure 9 because the required scaling 
factor to superpose the spectra is larger in this case since the accretion disk is not truncated. 
Note that there is a large discrepancy from the FUSEl spectrum. The 12,000K synthetic 
spectrum provides the best fit to the HST and FUSEl spectra, but the corresponding fit 
in the optical region is poor. 
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Fig. 16. — Models with four different values of WD T^s and an accretion disk that is 
isothermal in the inner region and follows a standard model in the outer region. The accretion 
disk temperature profile is in Table 5. The models fit the discovery SDSSJ0809 spectrum. 
The other observed spectrum is from the archive. The models are insensitive to the WD Teff 
at optical wavelengths and are essentially indistinguishable in the figure. The 45,000K WD 
contribution is at the bottom. 
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Fig. 17. — As in Figure 16 except in the UV. From top to bottom, the synthetic spectra 
show the models containing a WD with %s = 50,000K, 45,000K, 40,000K, and 35,000K. 
The lowest synthetic spectrum is the contribution of the accretion disk only. 
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Fig. 18. — Models with four different values of WD T^s and an accretion disk whose 
temperature profile is in Table 6. The models fit the SDSSJ0809 spectrum from the archive. 
Compare with Figure 16. The lowest curve is the contribution of a 45,000K WD. 
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Fig. 19. — As in Figure 18 but showing the UV fit. From top to bottom, the curves are 
models containing a WD with Teg = 50,000K, 45,000K, 40,000K, and 35,000K. See the text 
for a discussion. The lowest synthetic spectrum is the accretion disk contribution only. 



